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This  feasibility  study  concluded  that  a  computer  simulation  of  ship 
underway  replenishment  can  be  developed,  A  system  of  linear  equations 
of  motion  for  describing  two  ships  maneuvering  on  essentially  parallel 
courses  was  chosen  as  the  basic  equations  for  the  mathematical  model 
to  be  used  in  the  computer  simulation.  The  derivation  of  these  linear 
equations  1b  presented.  A  better  representation  using  nonlinear 
equations  of  motion  will  be  considered  in  a  later  phase  of  the  study. 
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DEFINITIONS  AND  TERMINOLOGY 


o,  x,  y,  z 


V  V  yr  V 


2'  2'  y2 ' 


Coordinate  system  which  moves  with  the  constant 
average  velocity  of  the  two  ships. 

Coordinate  axes  fixed  in  ship  1. 

Coordinate  axes  fixed  in  ship  2. 


u,  v 


-  Velocities  along  x  and  y,  respectively. 

-  Angle  of  yaw  measured  between  coordinate  axes 
ox  and  ox . 


M 


-  Mass  of  the  body  (ship) . 


I 

z 


-  Moment  of  inertia  of  the  ship  about  coordinate 
axis  oz. 


X,  Y 

N 


-  Force  acting  on  the  ship  along  coordinate  axes 
ox  and  oy. 

-  Moment  acting  on  ship  about  coordinate  axis  oz. 
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INTRODUCTION 


The  operational  procedure  of  replenishing  ships  at  sea  by 
steaming  on  parallel  courses  in  close  proximity!'^  cair.a  into 
general  use  during  World  War  II.  This  procedure  is  still  t  id 
extensively  by  the  Navy,  although  collisions  have  occurred  during 
repieniahment  operations  which  have  resulted  in  d, image. 

Analysis  of  the  replenishment  operations  may  provide  a  means 
of  reducing  the  collision  hazard  by  improved  methods  of  ship 
maneuvering  and  ship  control. 

The  dynamic  interaction  of  two  ships  in  close  proximity 
maneuvering  on  parallel  courses^-  involves  the  hydrodynamic 
interaction  between  the  ships,  the  skills  of  the  helmsmen,  and 
the  conning  officer  of  both  ships.  This  complex  dynamic 
interaction  is  not  completely  understood  and  requires  further 
analytical  and  experimental  study. 

This  laboratory  is  currently  engaged  in  an  applied  researcn 
program  to  define  the  parameters  affecting  ship  control  during 
underway  replenishment  (UNREP) .  Monitoring  of  these  control 
parameters  aboard  ship  may  provide  information  for  the  conning 
officer  and/or  helmsmen  which  will  reduce  the  collision  hazard 
and  also  improve  the  efficiency  of  the  UNREP  operation. 

This  first  phase  report  from  the  UNREP  research  program 
describes  the  background  of  the  problem,  derives  the  linearized 
equations  of  motion  for  two  ships  on  essentially  parallel  courses, 
and  discusses  the  computer  simulation  to  be  developed  at  NSRDC 
from  these  equations  as  uhe  basis  of  the  mathematical  model. 

Future  work  in  the  UNREP  program  will  concentrate  on  the 
development  of  the  computer  simulation  using  the  linear  equations 
of  motion,  analysis  of  the  computer  simulation  results,  and 
refinement  of  the  mathematical  model  to  include  nonlinear  terms. 

It  may  be  necessary  to  go  directly  to  the  nonlinear  maneuvering 
equations  if  preliminary  investigations  indicate  the  presence  of' 
instabilities  in  the  solutions  of  the  linear  equations. 


^Superscripts  ■'ter  to  similarly  numbered  entries  in  the  Technical 
References  at  chc  end  of  the  text. 
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GENERAL  DISCUSSION  OF  UNDERWAY  REPLENISHMENT 


Replenishment  operations'^'^  are  conducted  at  sea  for  the 
purpose  of.  transferring  stores  between  ships  to  extend  the  period 
of  time  that  the  ships  can  be  operational  at  sea.  Since  the 
cargo  must  be  guided  and  controlled  during  the  transfer  operation, 
a  suitable  physical  connection  must  be  maintained  at  all  times 
between  the  ships  moving  essentially  parallel  to  each  other  and 
at  nearly  identical  speeds.  This  physical  connection  requires 
that  the  ships  be  close  together;  thus,  there  is  the  danger  of 
collision.  The  tracking  ship  is  usually  assigned  the  task  of 
avoiding  collision  and  maintaining  station  relative  to  the  leading 
ship.  The  leading  ship  must  maintain  a  steady  course  and  must 
keep  oscillations  about  this  course  to  a  minimum. 

During  UNREP  operations  both  relative  speed  and  separation 
distance  of  ships  are  monitored  by  the  conning  officer4-  on  the 
tracking  ship.  A  line  with  markers  is  used  to  measure  the 
distance  between  the  ships.  One  end  of  the  line  is  secured  to 
the  leading  ship;  the  other  end  is  te.'.ded  by  men  on  the  tracking 
ship  who  take  up  and  pay  out  the  line  to  keep  it  taut.  The 
conning  officer  orders  small  course  and  small  speed  changes  to 
maintain  position.  Figure  1  shows  a  control  block  diagram  of  the 
UNREP . operation  with  the  tracking  ship  on  the  starboard  beam  of 
the  leading  ship.  This  diagram  is  discussed  in  appendix  A. 


Deiire* 

count 


2 

From  Thal-T. arsen. 


Figure  1  -  Replenishment- At  Sea  Operation 
with  Tracking  Ship  on  Starboard  Hearn 
of  Leading  Ship 
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The  constant  tensioned  wire  highlines  installed  in  the 
Navy's  newly  constructed  UNREP  ships  and  in  UNREP  ship  conversions 
introduce  another  parameter  to  be  considered.  This  parameter  is 
a  steady  lateral  pull  due  to  constant  tensioned  wire  highlines 
between  the  two  ships  which  may  affect  the  steering  of  both  the 
leading  ship  (delivery  ship)  and  tracking  ship  (receiving  ship) . 
However,  in  UNREP  operations-*  the  constant  tensioned  wire  high- 
lines  did  not  introduce  a  significant  problem.  This  area  is 
discussed  in  appendix  B. 

REPLENISHMENT  SIMULATION  TO  BE  DEVELOPED 

Pertinent  literature  relating  to  UNREP  opei  ions  was 
reviewed  to  evaluate  the  feasibility  of  computer  simulation  of 
underway  replenishment.  The  selected  papers  and  reports  from 
this  literature  search  were  studied  from  two  aspects:  their 
applicability  to  a  good  mathematical  model  of  UNREP  operation 
(e.g.,  assumptions  made,  the  reasons  for  these  assumptions,  and 
the  purpose  of  the  papers  and  reports) ;  and  their  applicability 
from  the  standpoint  of  simulating  ship  replenishment  by  computer 
(e.g.,  selection  of  the  computer,  the  reasons  for  its  selection, 
its  capabilities  and  deficiencies,  and  its  availability). 

ANALYSIS  OF  PAST  WORK 

In  the  field  of  UNREP,  the  literature  is  scarce.  Newton,** 
Silverstein , 5  and  Taylor6  report  experimental  and  theoretical 
investigations  of  interaction  forces  between  two  ships  maneuvering 
on  parallel  courses.  The  control  problem  involving  the  conning 
officer  and  helmsmen  during  an  UNREP  operation  is  not  considered 
in  these  investigations. 

Methods  presented  in  these  references  m^y  be  helpful  in 
determining  the  forces  between  the  two  ships  during  the  UNREP 
computer  simulation  in  this  program. 

The  University  of  Calif ornia ^ ^  reported  significant  UNREP 
work  clone  10  years  ago  by  a  research  group  sponsored  by  NAVSHIPS. 
Thic  group  recorded  incomplete  statistical  data  during  full-scale 
ship  UNFEP  operations.  The  sea  trials  consisted  of  refueling  a 
destroyer,  USS  HALSEY-POWELL  (DD  686),  by  the  USS  ASHTABULA 
(AO  51)  at  17  knots,  and  replenishing  USS  UHLMANN  (DD  687)  by  the 
USS  TICONDEROGA  (CVA  14)  at  25  knots.  The  project  was  terminated 
by  the  Navy,  however,  before  all  data  for  a  complete  statistical 
analysis  could  be  obtained.  Details  of  these  sea  trials  and 
exanples  of  these  data  appear  in  appendix  C. 
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Figure  1  shows  the  control  block  diagram  that  was  designed 
for  the  partially  competed  UNREP  analog  computer  simulation  at 
the  University  of  California.  This  block  diagram  represents  the 
UNREP  operation  with  the  tracking  ship  on  the  starboard  beam  of 
the  leading  ship.  The  incomplete  sea  trial  UHREF  data  mentioned 
above  were  uaed  to  determine  the  transfer  functions  incorporated 
in  tho  block  diagram.  A  brief  discussion  of  the  determination  of 
the  transfer  function?  appears  in  appendix  A. 

This  University  of  California  work  provideo  important  back¬ 
ground  information  useful  in  understanding  the  relationship 
between  ship  control  and  ship  interaction  during  UNREP  operations. 
However,  the  method  of  UNREP  computer  simulation  considered  by  tho 
University  ofi  California  was  not  suitable  for  use  in  this  program 
for  tho  following  reasons? 

»  Funds  in  excess  of  tho  amount  presently  allocated 
would  be  required  to  obtain  complete  statistical  data  for  this 
simulation  from  extensive  UNREP  operations. 

•  If  models  tests  are  used  to  obtain  UN REP  data,  the 
scaling  problems  involved  with  human  operators  in  the  control 
loops  for  those  tests  are  extremely  difficult  to  handle. 

•  This  analog  computer  simulation  using  linear  transfer 
functions  is  uxtrumeiy  hard  to  generalise  to  the  nonlinear  case. 

'ailing  and  Wood7,  present  the  linear  equations  of  motion 
for  two  ships  maneuvering  on  essentially  parallel  coursos.  It 
was  decided  to  une  these  equations  as  a  basis  for  the  mathematical 
model  in  tho  unrep  computer  simulation  to  bo  developed  for  tho 
following  reasons i 

•  Costly  hydrodynamic  work  con  he  eliminated  because 
most  of  tho  hydrodynamic  coefficients  in  tho  equations  of  motion 
can  bu  estimated  from  theoretical  considerations. 

e  Nonlinear  terms  can  be  addud  to  the  equations  of 
motion  in  a  later  phase  ol'  the  work  if  it  is  found  necessary  to 
consider  nonlinear  aspect a  of  UNREP  operations. 

LINEAR  EQUATIONS  OF  MOTION 

The  following  assumptions  are  made  in  developing  the  linear 
equations '  of  motion. 
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•  The  ships  operate  on  nominal  parallel  courses;  thus, 
the  time  average  resultant  velocity  vectors  are  parallel. 

e  The  average  speed  of  both  ships  is  constant  and  equal. 

e  The  short-time  deviations  from  the  above  conditions, 
are  small.  Thus,  the  linear  displacement  from  the  average 
posicion  of  the  ship  is  small  compared  to  the  ship  dimensions; 
angular  deviations  are  small;  velocity  disturbances  are  small 
compared  to  the  average  velocity;  and  accelerations  corresponding 
to  these  disturbances  result  in  moments  and  forces  which  are 
small  compared  to  those  already  acting  on  the  ship. 

In  order  to  describe  the  motion  of  the  two  ships  it  io 
necessary  to  use  two  sets  of  coordinate  systems, ^  shown  in  figure  2. 


From  Hauling  and  Wood. 

Figure  2 

Coordinate  Systems 
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One  coordinate  system  (o,  x,  y,  z)  moves  with  the  constant  average 
velocity  of  the  two  ships,  oz  is  directed  vertically  downward; 
the  xy  plane  coincides  with  the  surface  of  the  water;  and  ox  is 
parallel  with  the  average  velocity  of  the  two  ships.  The  second 
set  of  coordinates  (O^Xj/y^zi)  and  (o2  ,x2 ,V2  •  z2>  is  fixed  in 
the  two  ihips  with  o  fixed  at  the  ship's  center  of  gravity,  ox 
directed  forward,  oy  to  starboard,  and  oz  downward.  The  plane 
oxz  is  the  ship's  plane  of  longitudinal  symmetry.  The  plane  oxy 
is  parallel  to  the  water  surface  when  the  ship  floats  at  rest. 

Zt  is  assumed  that  the  principal  axes  of  inertia  are  ox,  oy,  and 
oz.  This  assumption  is  exactly  correct  for  oy  and  nearly  correct 
for  ox  and  oz  where  ships  of  normal  proportions  and  mass  distri¬ 
butions  are  considered.  The  motion  of  the  ship  is  restricted  to 
the  xy  plan«»;  and  thus,  pitch,  heave,  and  roll  motions  are  neglected. 
The  motion  is  described  by  velocities  ,  u2,  vj_ ,  vj,  acceler¬ 
ations  di,  u2,.<ri,  <12  parallel  to, the, x  and  y  axes,  and  angular, 
displacements  Vj.,  V2,  velocities  tj.,  T2,  and  accelerations  ^1,  V 2 
about  the  z  axis. 

In  the  two  coordinate  systems7 fixed,  respectively ,  in  the 
two  ships,  the  equations  of  motion  are  described  as  follows  (the 
subscripts  1  and  2  rofer  to  ships  one  and  two,  respectively) 1 
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where 


M  ®  mass  of  the  ship 

I  *  moment  of  inertia 

z 

T  a  thrust  of  propellers 

Xr,Yr  =>  rudder  forces  applied  to  the  ship 

Nr  *  rudder  moment 

X,Y  *  motion-dependent  external  hydrodynamic  forceB 

applied  to  the  ship 

N  «  motion-dependent,  external  hydrodynamic  moment 

applied  to  the  ship. 

It  is  assumed  that  the  hydrodynamic  variables7  (X. ,  X_,  Y  , 
*2'  Nl'  dopond  on  the  followings  1 


Velocities  u. .  u„,  v,  ,  v„ 

X  *  0 .  d, 

Accelerations  ii^,  Cs2,  , 


'  *1'  ^2  for  each  ship. 

s  I  It 

v2,  *2  for  each  ship. 


Relative  positions  for  each  ship 


xx  and  y2  -  y^ . 


Thus,  the  variable  X^  may  be  expressed  ao^X^  "_xl(ui»  uz •  VI»  v2 » 

u]/  ^2'  ^1'  V2'  *^1'  ^2'  ^1'  ^  2.'  x2  **  Ki  '  ^2  *  ’/l5  w  8‘mflflr 
expressions  for  X2,  YA,  Y2,  N.,  N_.  Other  external  forces  not 
mentioned  might  result  from  wind,  wavrs,  etc.  All  of  these  forces 
and  moments  are  dependent  on  the  geometry  of  the  ship. 

It  is  assumed  here  that  the  motion  of  each  ship  consists  of 
small  perturbations  superimposed  on  an  initial  state  of  motion.7 
Thus, the  perturbation  of  the  u  velocity,  for  example,  is 
represented  by 


u(t)  -  uQ(t)  *■  6u(t) 


(?) 
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where  6u  is  a  small  quantity.  Assuming  the  X2  is  an  analytic 
function  of  the  independent  variables,  X2  is  expanded  in  a 
Taylor's  series  about  its  value  at  the  initial  state  of  motion: 

Xi  =  Xl(ulo'  u20 . (y2  “  5l>o 


3X.  2X 

+  r —  6u.  +  r —  6u_,  +, 
BUj^  1  3u2  2 


3X, 


.  + 


3(y2-  yx) 


5 (y,-  y. ) 


(3) 


i  ’S  ? 

+  5  -T ~  <4ui>  + 

3  u. 


Substituting  equation  (2) ,  plus  other  equations  representing 
the  perturbations  of  its  motion,  and  equation  (3)  into  the  first 
of  equation  (1)  and  retaining  terms  to  the  first  power  of  6u^  etc 
results  in  the  equation  for  the  initial  state  of  motion,7 


Ml(u10"  V'lO1  -  Vu10'  U20'--'  <V  'l1*1  *  T10+  XR10 


plus  the  equation  for  perturbed  motion 


(4) 


M  (6u  +  Y  6v  -  i!  v  ) 

1V  1  10  1  1  lo' 


3X  3X.  3X 

5—  sUl  ,  j-  «u2  ♦  ...  — - — 

i  2  s  (y2-  y2) 


f  .y2-  yj)  +  «V  4xm 


where  tho  nonlinear  terms  are  dropped; 

■i>2)  • 

These  nonlinear  terms  will  be  considered  in  future  work. 


(5) 


/  3 

(  l  ‘I 

^example  ^  -y—  •  (5u, 
3  u^ 
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It  should  be  noted  that  the  partial  derivatives  3X^/'du^  etc 
are  evaluated7  at  the  initial  state  of  motion  (i!]_  *  uiO»  U2  3  U2Q, 
y2~  y,  =  (y2-  Vi^o1  The  changes  in  thrust  and  rudder  force,  <$T 
and  6x^  are  considered  to  depend  on  the  motion  or  to  be  independ¬ 
ent  variable  quantities.  However,  it  is  required  that  they  be 
email  quantities. 

It  is  aljo  necessary  to  introduce  the  assumgtion  of  constant 
average  velocity,  uQ,  of  both  ships  parallel  to  ox.  The  resultant 
velocities  ir.  the  cx  and  oy  directions  are  givon  in  terms  of  tne 
velocities  parallel  to  ox  and  oy  for  either  ship  as  follows: 


u  =  u  cos  ?  +  v  sin  f 

v  ^-u  Bin  V  +  v  cos  V  ,  (6) 


where  the  subscripts  denoting  either  ship  have  been  dropped. 

By  aubfltitutir  the  equations 

u  ■»  u  +  6u 
o 

V  ■  v  ■»  4  v 

o 

u  ■  u  +  5u 
o 

v  ■  Sv 

V  -  Y  +  6f  (7) 

c 


into  equation  (6) ,  the  following  two  equations 


7 


are  obtained: 


u  +  <$u  -  (u  +  6u)  cos  (Y  +  6Y)  ♦  6v  sin  +  6Y) 
000  o 


v  +  6v  «  -(u  +  iu’;  sin  (Y  +  <5Y)  +  4v  cos  <Y  +  6Y)  .  (B) 

000  o 
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It  is  assumed  here  that  the  average  yaw  (y)  angle  of  the  ship 
is  small7  so  that  the  following  approximations  may  be  made 


cos  y  h  l 
o 

Bin  y  2  ¥  .  (9) 

o  o 


Then  substituting  these  vproximations  (equation  (9))  into 
equation  (8)  results  in  the  equations, 


u  +  4u  =  u  +  (6u  +  y  6v) 
o  o  o 


v  +  5v  ■ 
o 


-y  u  -u  ay  +  6v 
o  o  o 


(10) 


Thus,  to  the  first  order  approximation 


7 


u  - 

o  o 


u 


v  ■  ~y  u 
o  o  o 


(11) 


where  uQ  is  the  constant  average  velocity  of  both  ships  and  f  0 
the  appropriate  average  angle  of  yaw  for  the  ship  in  question. 
Since  n'0  has  been  assumed  to  be  small,  terms  of  the  form  y0iu, 
etc,  can  bo  neglected. 

The  equations  of  perturbed  motion7 about  an  initial  steady- 
otato  motion  for  the  two  ships  are  (uq  ■  constant,  v^  »  constant)! 
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t 


3X  3X  3X 

HA  ■  5^iui *  ir2  6*2  *  ■■■  *  — 1 


3  (y2-  y:) 


<(V  *1>  +  *V  SXR1 


ay,  ay. 


W  \  V  ’  aTT  Sui  *  aST 


3Y, 


+  .  .  .  ■  + 


3 (y2-  Yj) 


4(y2-  yt)  +  «yR1 


3N 


1 


3N , 


‘il’i  ‘  5^  4ui  +  * 


a(y2-  y2i 


4  (y2~  yx) 


3x2  3X  3  X 

m2u2  "  FuT  4u:  +  ‘JuT  4  u2  +  *  •  *  77  ^7 

1  2  3(y2-  yx) 


4  (y2-  y2) 


+  6T2  +  iXR2 


•  _  3^  3y 

M,(v,  +  V,uJ  -  37i  4u,  +  — ±  3u_ 

2  2  2  o  3u^  1  2 


3Y. 


+  ...  — - - My,-  y. )  +  <5 Y 


3  (yj-y^ 


R2 
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J 


3N2  3N2  .  -  - 
1z2,2  ’  ST  4U1  +  •"  *  -r-r—  s<y2-  y  ) 

1  3 <y,-  y,) 


(12) 


The  retention  of  only  tne  first  order  hydrodynamic  force 
terms  of  equation  (12) ^  implies  that  these  forces  can  be  expressed 
by  terms  proportional  to  the  first  powers  of  perturbation 
displacements,  velocities,  and  accelerations.  However,  Wehausen8 
and  Cummins'^  shewed  that  the  hydrodynamic  forces  on  a  body  moving 
at  a  free  surface  can  be  expressed  by  the  first  order  perturbation 
terms  only  in  special  cases.  One  such  c<*.el°  is  that  of  small 
amplitude  periodic  motion,  for  example,  that  of  ship  motion  in 
regular  waves  where  the  amplitudes  of  the  motion  are  small.  The 
present  problem  will  not  necessarily  be  confined  to  the  case  of 
periodic  motion.  Thus,  it  is  uncertain  that  the  retention  of 
only  linear  terms  sufficiently  describes  the  hydrodynamic  forces 
involved . 


Accordingly  the  computer  simulation  using  the  linear 
maneuvering  equations  may  exhibit  unstable  solutions  thus  making 
it  necessary  to  go  directly  to  the  nonlinear  maneuvering  equations. 
The  nonlinear  equations  will  be  obtained  by  adding  the  perturbed 
higher  order  force  terms  to  the  linear  equations. 


COMPUTER  FACILITY 


The  computer  simulation  for  the  UNREP  project  must  allow 
for  a  man-in-the-ioop.  To  obtain  this  man-machine  interface  the 
simulation  must  operate  in  real  time,  thus  implying  that  the 
simulation  be  done  on  an  analog  computer.  This  laboratory 
presently  has  an  analog  computing  capability  which  should  be 
sufficient  for  the  simulation.  This  capability  includes  the 
interface  to  a  mock-up  of  a  surface  ship  bridge.  This  facility 
is  being  enlarged  to  include  a  complete  hybrid  computing 
capability  which  will  allow  the  simulation  to  run  in  real-time 
while  providing  the  programing  flexibility  of  a  digital  computer. 

DISCUSSION 

There  is  very  little  quantitative  data  on  the  two-ship  hydro- 
dynamic  coefficients  available;  therefore,  it  will  be  necessary 
to  estimate  these  hydrodynamic  coefficients  for  a  general  ship 
situation  in  the  maneuvering  equations.  The  UBe  of  estimated 
coefficients  is  justified  to  some  extent  because  the  study  is  not 
directed  toward  specific  ships  but  toward  general  results 
applicable  to  a  variety  of  ship  types.  These  coefficients  will 
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be  estimated  for  the  most  part  from  unpublished  work  from  the 
hydrodynamics  facility  of  NSRDC.  Every  effort  will  be  made  to 
use  hydrodynamic  coefficients  which  specifically  refer  to  the 
two-ship  situation.  When  these  coefficients  are  not  available, 
the  single-ship  hydrodynamic  coefficients  will  be  used. 

The  first  studies  made  by  exercising  the  UNREP  simulation 
will  attempt  to  determine  the  minimum  number  of  maneuvering 
parameters  necessary  for  use  by  the  conning  officers  and/or 
helmsmen  for  station  keeping  during  steady-state  UNREP  operations. 
The  maneuvering  parameters  for  steady-state  operation  may  not  be 
sufficient  for  accurate  station  keeping  during  non-steady-state 
operations. 

These  first  studies  will  also  include  determining  the  display 
characteristics  of  the  maneuvering  parameters  and  investigating 
the  accuracies  needed  for  the  sensors  measuring  these  parameters 
aboard  ship. 

After  implementing  the  UNREP  computer  program  or  the  hybrid 
computer  the  results  will  be  compared  with  the  limited  available 
unpublished  and  published  'JNREP  data  from  sea  trials  and  model 
tests.  This  validation  procedure  will  also  include  parametric 
studies  of  the  general  hydrodynamic  sensitivities.  As  a  result 
of  these  sensitivity  studies,  it  may  bo  found  necessary  to 
determine  some  hydrodynamic  coefficients  more  accurately  by  model 
tests  at  NSRDC. 

Other  areas  of  future  investigation  include; 

s  The  adaption  of  the  UNREP  compucor  simulation  to  the 
ship  approach  situation  at  the  beginning  of  the  UNREP  maneuvers . 

•  Modification  of  the  UNREP  computer  simulation  to 
ir  'de  sea  3tate,  roll,  pitch,  etc. 

r  Collision  course  avoidance  and  the  display  of 
maneuvering  par  .'.meters  to  warn  the  conning  officer  and/or  helmsmen 
of  the  collision  situation  and  provide  information  to  reduce  the 
hazard. 


CONCLUSIONS 

From  the  limited  references'-^  available  on  underway 
replenishment,  it  was  determined  that  the  best  basis  for  a  mathe¬ 
matical  model  for  the  simulation  of  underway  replenishment  was 

the  line  equations  developed  by  Pauling  and  Wood  which  can  be 
general,  'cd  to  the  nonlinear  case. 
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Some  of  the  general  ship  hydrodynamic  coefficients  in  the 
equations  of  motion  may  have  to  be  determined  by  model  tests. 

The  method  to  simulate  underway  replenishment  using  classical 
control  theory  as  partially  developed  by  the  University  ^f 
California  is  not  suited  to  this  work  because  of  the  expense  in 
obtaining  statistical  data  used  in  the  simulation  and  because  the 
nonlinear  aspects  of  the  control  problem  are  extremely  difficult 
to  include. 

The  computer  program  to  be  developed  will  provide  important 
control  parameters  for  display  aboard  ship.  These  displays  will 
aid  the  conning  officer  and/or  helmsmen  in  ship  control  and  should 
reduce  the  risk  of  collisions  during  underway  replenishment 
operations . 
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APPENDIX  A 

STATISTICAL  ANATYC'  "  OF  SHIP  REPLENISHMENT  DATA 
REFERENCES 

(a)  Tichvinsky,  L.  M. ,  and  H.  Thal-Larsen,  "Replenishment  at 
Sea,  Experiments  and  Comments,"  1st  tech  rept  issued  under 
the  project  "Replenishment  at  Sea,"  Univ.  of  Cal? fornia, 
Berkeley,  Inst,  of  Engr.  Research  (now  Office  of  Research 
Services),  Rept  154,  Issue  1  (1  June  1960) 

(b)  Tnal-Larsen,  H. ,  "Manual  Steeriny  of  Two  Ships  in  Close 
Proximity  During  Replenishment-at-Sea  Exercises,"  contributed 
by  the  Automatic  Control  Div.  of  ASME  for  presentaion  at 
ASME  winter  annual  meeting,  Los  Angeles,  Cal.  (16-20  Nov  1969) 

This  appendix  presents  a  review  of  the  analysis  ~f  the  ship 
replenishment  data,  references  (a)  and  (b) ,  summarized  in  appendix 
C  and  carried  out  by  the  University  of  California. 

The  block  diagram  (figure  1  of  the  text)  shows  the  six 
variables  that  are  of  most  interest.  These  variables  on  the 
leading  ship  are  defined  as  o  =  rudder  angle  and  8^  =  deviation 
of  gyro  angle  from  the  desired  course.  On  the  tracking  ship,  6  = 
rudder  angle,  02  “  deviation  of  gyro  angle  from  the  desired  value, 
x  »  deviation  from  the  desired  separation  distance  of  the  two 
ships,  and  c<.  a  the  ordered  change  in  course  of  the  tracking  ship. 

In  addition  to  these  variables,  also  recorded  during  the  full- 
scale  ship  trials  were  the  helm  angles  on  both  the  leading  and 
tracking  ships. 

The  following  sign  conventions  were  established  for  the 
variables  (figure  1-A) . 


From  Tichvinsky  and  Thal-Larsen  (reference  (a)). 
Figure  1-A  -  Two  Ships  Abeam 


/fa 
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Since  compass  readings  increase  in  magnitude  when  a  ship  turns  to 
the  right,  the  gyro-angle  variables  0^  and  02  are  considered  co 
be  positive  for  a  ship  turning  to  the  right.  The  rudder-angle 
variables  a  and  6  are  assigned  positive  values  for  the  right- 
rudder  deflections.  An  increase  in  separation  between  the  two 
ships,  x,  is  taken  as  being  in  the  positive  direction.  The 
conning  officer  and  helmsmen  produce  negative  feedback  corrections. 
Thus,  a  will  be  negative  when  x  is  positive. 

Analysis  of  the  data  (appendix  C,  figures  1-C  and  2-C  was 
carried  out  to  determine  the  linear  transfer  functions  (figure  1 
of  the  text)  . 

The  17-knot  data  (appendix  C,  figure  1-C)  were  analyzed 
first.  Figure  1-C  shows  the  approximate  phase  relationship 
between  the  rudder  and  gyro  oscillations  for  the  oiler  (AO) , 

USS  ASHTABULA,  the  leading  ship.  The  gyro  was  approximately  160 
degrees  out  of  phase  with  the  rudder.  It  was  determined  that  ^he 
transfer  function  ip  nearly  equal  to 


_ 4 

D2(0.5D  +  1) 


(gyro  degree/min 


2 


* 

) /rudder  degree, 


(A-i) 


where  D  *  d/dt  is  the  differential  operator. 

Square-wave- rudder  tests,  at  different  frequencies,  applied 
to  the  ship  operating  alone,  confirmed  •  e  K]/D2  form  of  equation 
(A-l) .  The  additional  time  constant  in  equation  (A-l)  is  possibly 
caused  by  the  destroyer  on  the  starboard  beam  (references  (a)  and 
(b))  . 


For  the  helmsman  on  the  oiler  (AO)  it  iB  necessary  t_  find 
a  transfer  function  relating  the  AO  gyro  curve  as  input  to  the  AO 
rudder  curve  as  output.  The  inverse  of  the  previously  determined 
transfer  function,  equation  (A-l) ,  would  q  .ite  accurately  satisfy 
the  mathematical  relationship;  .he  resulting  third  order 
derivative,  however,  would  not  be  physically  descriptive  cf  the 
helmsman . 


•Abbreviations  used  in  this  text  are  from  the  GPC  Style  Manual, 
1967,  unless  otherwise  noted. 
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A  statistical  description  of  the  helmsman  is  provided  by  the 
correlation  functions  in  figure  2-A. 


From  Thai -Larsen  (reference  (b)). 


Figure  2-A 

Statistical  Input-Output  Relationship 
for  the  Helmsman  plus  Rudder  Servo 
on  the  Leading  Ship  (AO  51)  at  17  Knots 


The  smoothing  action  of  the  statistical  prococa  involved  in 
calculating  these  curves  produced  simple  curvos.  The  method  by 
which  these  cnrvos  were  computed  and  their  interpretation  as 
equivalent  input-output  functions  is  discussed  In  roforenco  (b)  , 
It  io  sufficient  to  state  here  that  the  solid  curves  in  figure 
2-A  can  bo  considered  as  the  input  an  lo  to  the  helmsman  and  the 
dashed  curve  the  helmsman's  response. 

From  figure  2-A,  it  is  noted  that  when  the  ship  swings  fro.? 
maximum  velocity  from  one  side  to  the  other,  the  helmsman  applied 
maximum  rudder  corrr ition  to  stop  that  swing.  Thi i  type  of 
response  is  termed  i.  -rivativu  response.  Actual  1*,  the  helmsman  is 
about  10  seconds  1 a t •  j ,  which  is  about  the  time  . t  takes  h  m  by 
watching  the  gyro  repeater  to  estimate  when  macimum  velocity  has 
been  reached,  Thus,  the  proposed  transfer  function  model  for  tl;a 
helmsman  and  ruddar  servo  combination  becomes 


K 


2 


De"LlD 


rudder  degree/gyro  degrer 


( A — 2 ) 


!% 
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whore 


•  1  minute 


"1/6  mlnuto  or  10  seconds. 

Tho  mlnuo  sign  accounts  for  the  rudder  action  being  applied  to 
oppose  the  swing  of  the  ship. 

The  statistical  curves  in  figure  3-A  may  ba  regarded  as 
input  and  output  curves  for  the  conning  officer  on  the  tracking 
ship'.. 


Iron  Thai - h«r» en  (nfir«nu  (6)). 


Figure  3-A 

fltatisUcal  Input-Output  Relationship 
for  the  Conning  Officer  on  the 
Tracking  I'hip  (DD  086)  at  17  Unote 


it  was  determined  thnt  his  transfer  function  is  approximately 


-K 


D 


degree  course  olwinge/f  net  sopnration  distance  change 


( A — 3 ) 


•  /V 
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K.,  ■  0.04  degree/feet 


L2  -  1/12  minute  or  5  seconds. 

K,  was  estimated  by  taking  the  ratio  of  the  two  central  peaks  of 
the  curves  in  figure  3-A. 

Unfortunately,  Incomplete  data  allowed  only  rough  estimates 
to  be  made  of  the  transfer  functions  for  the  conning  officer  and 
the  helmsman  of  the  tracking  ship.  Tor  simplicity  it  was  assumed 
that  the  helmsman  on  the  DD,  together  with  the  rudder  servo,  wore 
represented  by 


rudder  degree/gyro  degree 


(A-4 ) 


where  K4  equals  approximately  two  in  the  frequency  range  of 
Interest  displayed  in  figure  1-C  of  appendix  C.  Evidence  from 
reference  (a)  and  from  laboratory  tests  (reference  (b) )  indicated 
a  transfer  function  for  the  tracking  ship  (DD)  of  the  formt 


hre~L3D 

5 

L (TD  +  1) 


(gyro  degree/min) /rudder  dogroo 


(A-5) 


where 


T  £  1/12  minute  or  5  seconds 
L3“  1/15  minute  or  4  seconds 

Krt  6  minutes"1 
5 

The  2b-knot  data  in  figure  2-c  were  analysed,  it  is  evident 
that  the  holmsman  operates  the  rudder  at  a  more  rapid  rate 
compared  to  the  17-knot  case.  These  statistical  curves  are  used 
for  determining  the  char aotoristic  properties  of  the  dynamic 
system, 

1 
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Figure  4-A  shows  the  atatiatical  input-output  curvea  for  the 
leading  ahip,  an  aircraft  carrier  (CVA) . 


Figure  4-A 

Statistical  Input-Output  Relationship 
for  the  Leading  Ship  (CVA  14) 
at  25  Knots 


The  two  curvea  are  approximately  180  degroee  out  of  phane.  It 
waa  aaaumod  that  a  linear  tranafar  function  ia 

K6 

gyro  degree/rudder  degree  (A-6) 

D 


2 

where  ia  approximately  0,14  (yyro  degree/minute  ) /rudder  degree. 

'’‘he  helmsman  doea  not  have  much  time  to  eatimate  angular 
velocitiea  of  the  gyro  at  the  higher  speed,  It  waa  aaaumod  that 
tho  transfer  function  of  the  helmsman  together  with  the  rudder 
servo  to  be  approximately 

rudder  degree/gyro  degren  (A-7) 


3ti 
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where  h7  varies  between  1.'  and  7.  Thus,  it  is  possible  for  the 
helmsman  to  compensate  for  the  sensitivity  of  the  ship.  This 
sensitivity  decreases  with  the  square  of  the  frequency. 

Figure  5-A  sr.^/s  the  self  spectral  density  curves  for  the 
rudder  and  the  gyro  on  the  leading  ship  (CVA  14)  at  25  knots. 


From  Thsl-Larsen  (reference  (b)). 
Figure  5-A 

Self -Spectral  Density  Curves  for  the  Rudder 
and  the  Gyro  on  the  Leading  Ship  (CVA  14) 
at  25  Knots 


This  figure  shows  that  the  aircraft  carrier  is  insensitive  to 
relatively  nigh-f roquency  rudder  perturbations.  While  the 
spectral  density  curve  for  the  rudder  shows  a  number  of  peaks 
at  frequencies  above  w  ■  1  radian/min,  no  such  peaks  are  seen  in 
the  epectral  density  curve  for  the  gyro  in  figure  5-A. 

The  source  of  the  strong  spectral  density  peaks  of  the  CVA 
rudder  at  circular  frequencies  of  about  10  and  15  radians/min 
(corresponding  to  periods  of  approximately  40  and  26  seconds, 
respectively)  is  not  clear.  These  frequencies  *re  evident  in  the 
solid  curve  in  figure  4-a.  Also  the  rudder  spe^  al  density 
curve  has  sm-.1 1  peaks  corresponding  to  frequencies  of  15  and  13 
seconds,  respectively.  The  40-  and  26 -second  periodo  are  possibly 
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induced  by  the  action  of  the  aea  and  the  wind.  The  15-  and  13- 
second  periods  are  possibly  due  to  the  rapid  and  continuous 
swinging  back  and  forth  of  the  rudder  which  gives  the  helmsman  a 
"feel"  for  the  ship. 

The  characteristics  of  the  helmsman  on  the  CVA  are  provided 
by  figure  6-A, 


From  Tha1-Larapn  (reference  (b)). 


Figure  6-A 

Statiotical  Input-Output  Relationship  for  the 
Helmsman  plus  Rudder  Servo  on  the  Leading 
Ship  (CVA  14)  at  25  Knots 


The  solid  curve,  representing  the  gyro,  does  not  extend  enough  to 
display  the  6-minute  period  of  continuous  oscillation.  If  the 
solid  curve  is  recorded  as  input  to  the  helmsman,  the  helmsman's 
response  has  a  gain  factor  of  two.  This  is  more  than  the  1.4 
from  equation  (A-7)  and  represents  a  deviation  in  the  correct 
dire-  .r-n  because  of  the  higher  frequency  involved  in  figure  6-A. 


JU 
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It  wao  assumed  that  the  trnnpfer  function  of  the  tracking 
ship  (DO)  would  probably  be  in  the  form 


-r  gyro  degree/rudder  degree  (A-B) 

D 

where  Kg  is  probably  less  than  200  (see  figure  7-A) . 


From  Thal-Laraen  (reference  (b)), 


Figure  7-A 

Statistical  Curve  for  the  Ruddx -  of  the 
Tracking  Ship  (DO  687)  at  25  'ots 


Lack  of  ship  trial  test  data  prevented  a  more  accurate 
estimate  of  equation  (A-8)  as  it  does  for  equation  (A-9)  for  the 
helmsman  on  the  tracking  ship. 


rudder  degree/gyro  degrea. 


(A-9) 
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It  was  determined  from  figure  8 -A  that  the  transfer  function 
for  the  conning  officer  is  essentially 


degree  course  change/feet  separation  distance 

(A-10) 

which  is  the  same  form  as  equation  (A-3)  for  the  17-knot  test. 

Kip  was  found  to  be  approximately  0.01  degree/feet  and  L4  nearly 
equal  to  1/12  minute. 


From  Thal-Larsen  (reference  (b)). 


Figure  8-A 

Statistical  Input-Output  Relationship  for  the 
Conning  Officer  on  the  Tracking  Ship  (DD  687) 

at  25  KnotB 


As  a  result  of  this  project  being  discontinued  by  the  Navy, 
only  incomplete  sea  trial  data  wore  available  for  this  analysis. 
Nevertheless,  the  statistically  reduced  data  obtained  may  be  of 
use  to  investigators  in  future  work. 
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APPENDIX  B 


MODERN  UNDERWAY  REPLENISHMENT  OPERATIONS 


This  appendix  gives  a  summary  of  the  study  made  by  Captain 
Lienhard  (USN) ,  commanding  Officer  of  the  USS  SYLVANIA  ( AFS  2)* 
on  the  effects  of  constant  tensioned  wire  highlines  used  during 
modern  underway  replenishment  operations. 

When  transferring  operations  are  conducted  by  conventional 
wire  highline,  housefall,  or  Burtoners,  there  is  a  transient 
effect  on  steering  as  the  load  passes  between  the  two  ships. 
Normally  loads  transferred  by  this  method  do  not  exceed  3,500 
pounds.  The  maximum  line  pull  will  not  greatly  exceed  this  amount 
if  a  proper  catenary  is  maintained.  However,  in  constant 
tensioned  transfer,  the  ram  tensioner  may  be  set  to  apply  the 
maximum  line  pull  of  12,500  pounds  in  the  stores  ship  (AFS), 

14,000  pounds  in  the  oiler  (AO),  or  18,000  pounds  in  the  ammu¬ 
nition  ship  (AE) ,  for  cargo,  fueling,  or  missile  transfer  station 
rigs  respectively. 

The  Commanding  Officer  of  the  USS  SYLVANIA  (AFS  2)  conducted 
a  survey  of  the  effects  of  constant  tensioned  lines  on  steering 
over  an  18-month  period  on  a  variety  of  different  types  of 
customer  ships.  The  Commanding  Officer  of  each  customer  ship  was 
then  asked  the  effects  on  his  ship’s  steering.  The  results  are 
summarized  as  follows: 

•  Use  of  constant-tensioned  wire  highline  produces  no 
greater  problem  to  safe  ship  handling  than  conventional  wire 
highline,  housefall,  or  Burtoning. 

•  Possioly,  the  hazard  may  be  less  for  the  constant- 
tensioned  wire  highline  because  the  tensioned  pull  is  a  constant 
parameter . 

In  the  opinion  of  the  author  the  principle  elements  of  tnis 
new  parameter  which  may  cause  a  steering  problem  for  the  leading 
(delivering)  and  tracking  (receiving)  ship  are: 


*Lienhard,  B.  A.  (Capt,  USN) .  "Shiphandling  in  Modern  Underway 
Replenishment,"  NAVSHIPS  Tech  News,  Vol.  21,  Ho.  4 
NAVSHIPS  01C0  000  2085  (Apr  1972) 


•  Amount  of  pulled  tension. 

•  Location  of  the  attachment  points  of  the  rigs  in 
rexation  to  the  location  of  the  ship's  pivotal  points. 

•  Draft  and  trim  of  the  ships. 

There  were  also  important  safety  considerations  which  were 
mentioned  by  Captain  Lienhard: 

•  The  tracking  ship's  conning  officer  can  hasten  the 
initial  hookup  and  improve  maneuvering  safety  by  coming  alongside 
of  the  leading  ship  at  no  greater  than  normal  replenishment  ship 
separation  distance  when  the  lines  and  rigs  are  being  passed 
between  the  ships.  Then  the  conning  officer  should  move  the 
tracking  ship  out  to  a  greater  separation  as  the  tension  is 
applied  in  the  correct  amount.  The  tracking  ship  can  adjust  the 
lateral  separation  to  suit  the  present  situation  of  speed,  wind, 
and  sea,  and  amount  of  rudder  being  carried. 

•  Once  tenaion  is  applied  on  the  lines,  the  conning 
officers  of  both  ships  should  be  alert  to  adviBe  each  other  imme¬ 
diately  of  a  developing  situation  which  might  require  an  emergency 
breakway.  For  personnel  safety  and  elimination  of  potential  ship 
equipment  damage,  the  leading  3hip  should  detension  the  rig  before 
t.ie  tracking  ship  attempts  to  release  the  pelican  hook  or  unshackle 
and  cast  off  the  rig. 
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APPENDIX  C 


FULL-SCALE  UNDERWAY  REPLENISHMENT  SEA  TRIALS 
REFERENCES 

(a)  Tichvinsky,  L.  M. ,  and  H.  Thal-Larsen,  "Replenishment  at  Sea, 
Experiments  and  Comments,"  1st  tech  rept  issued  under  the 
project  "Replenishment  at  Sea,"  Univ.  of  California,  Berkeley, 
Inst,  of  Engrg,  Research  {now  Office  of  Research  Services) , 
Rept  154,  Issue  1  (1  June  1960) 

(b)  Thal-Larsen,  H. ,  "Manual  Steering  of  two  Ships  in  Close 
Proximity  During  Replenish  ent-at-Sea  Exercises,"  contributed 
by  the  Automatic  Control  Div.  of  ASME  for  presentation  at 
ASME  winter  annual  meeting,  Los  Angeles,  Cal.  (16-20)  Nov 
1969) 

The  University  of  California  (reference  (a)  and  (b) )  carried 
out  full-scale  sea  trials  on  underway  replenishment.  These 
studies  were  not  complete  because  the  contract  was  terminated 
before  complete  statistical  data  could  be  taken.  The  sea  trials 
consisted  of  refueling  a  destroyer  USS  HALSEY-POWELL  (DD  686)  by 
the  USS  ASHTABULA (AO  51)  at  17  knots  and  replenishing  USS  UHLMANN 
(DD  687)  by  the  USS  TICONDERGA  (CVA  14)  at  25  knots.  This 
appendix  summarizes  the  methods  used  for  studying  the  replenish¬ 
ment  operations  and  gives  examples  of  data  from  these  trials. 

During  the  trials  helm  angle,  rudder  angle,  and  gyro  compass 

heading  for  both  ships  were  recorded  as  a  function  of  the  time 

for  the  exercise  while  the  two  ships  were  abeam.  Course  commands 
by  the  conning  officer  of  the  tracking  ship  and  the  separation 
distance  of  the  ships  luring  tho  exercises  were  also  recorded 
(See  figureB  1-C  and  2-C.)  These  eight  variables  wore  UBed  for 
auto-  and  cross-correlation . 

Statistical  analysis  of  the  data  yielded  information  on  the 
dynamic  responsp  lags  at  different  speeds  and  of  factors  affecting 
human  behavior.  In  addition  to  the  variables  mentioned  above, 
the  following  were  found  to  be  significant! 

•  Automatic  pilot  data. 

•  Propeller  rpm  of  both  ships. 

•  Pitch  and  roll  angles  of  both  ships. 

•  Environmental  conditions!  weather,  sea  state, 

wind,  etc. 
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After  analysis  of  the  above  mentioned  data  from  the  replen¬ 
ishment  sea  trial3,  the  following  summarized  conclusions  were 
made : 


«  The  problem  of  primary  importance  is  station  keeping 
with  cargo  transfer  of  secondary  importance.  Improvements  in  the 
station-keeping  and  cargo  handling  technique  are  needed. 

•  The  dynamic  interaction  between  the  conning  officer 
and  the  helmsman  is  the  Drimary  reapon  for  variations  in  separation 
distance  between  the  two  ships. 

•  The  dynamic  interaction  between  tha  helmsman  and  the 
ship  and  the  hydrodynamic  interaction  between  the  two  ships  must 
both  receive  consideration. 

•  The  increase  of  sensitivity  with  speed  in  the  man- 
ah.ip  feedback  loop  results  in  the  rapid  increase  in  perturbation 
frequency  and  the  corresponding  decrease  in  amplitude  with 
increased  speeds. 

•  The  mean  rudder  angles  required  by  each  ship  results 
from  the  hydrodynamic  interaction  between  the  two  Bhips. 

•  The  hydrodynamic  interaction  between  two  ships  abeam 
produces  an  additional  l/.g  in  the  response  of  the  ship  to  the 
rudder  which  makes  steering  more  difficult. 

•  The  lack  of  knowledge  by  the  helmsman  of  the  hydro- 
dynamic  interaction  between  the  ships  result  in  large  separation 
distances  at  higher  replenishment  speeds. 

•  A  mathematical  description  on  a  statistical  basis  of 
the  dynamics  of  two  ships  and  the  helmsmen  of  the  tv:o  ships  and 
the  conning  officer  of  the  tracking  ship  involved  xn  ship  replen¬ 
ishment  can  bo  derived. 

The  following  recommendations  were  offered: 

•  A  means  should  be  developed  for  instantaneous 
indication  of  the  relative  positions  of  the  two  chips  at  close 
proximity.  This  would  enable  the  conning  officer  to  sec  at  once 
if  the  two  ships  are  on  convergent  or  divergent  courses. 

•  Further  study  of  the  station-keeping  problem. 

•  Experimental  and  theoretical  determination  of  the 
hydrodynamic  interaction  between  two  ships. 
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From  Thal-Laraen  (reference  (b)). 


Figure  2-c 

25-Knot  Test  Data  Abeam 
(Leading  Ship:  USS  TICONDEROGA  (CVA  14); 
Tracking  Ship:  USS  UHLMANN  (DD  687)) 
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•  Study  of  the  control  characteristics  at  high  speeds 
of  two  ships  during  replenishment. 

•  Investigate  the  use  of  ctrip-chart  recorders  for 
helping  the  helmsman  during  station  keeping. 
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From  Thal-Lcu  sen  (reference  (b)). 


Figure  1-C 

2 7 -Knot  Test  Data  Abeam 
(Leading  Snip:  USS  ASHTABULA  (AO  51'; 
Tracking  Ship:  USS  HALSEY-POWELL  (DD  686)) 
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